Bosonization out of equilibrium 
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We develop a bosonization technique for one-dimensional fermions out of equilibrium. The ap- 
proach is used to study a quantum wire attached to two electrodes with arbitrary energy distribu- 
tions. The non-equilibrium electron Green function is expressed in terms of functional determinants 
of a single-particle "counting" operator with a time-dependent scattering phase. The result reveals 
an intrinsic relation of dephasing and energy redistribution in Luttinger-liquids to "fractionalization" 
of electron-hole excitations in the tunneling process and at boundaries with leads. 

PACS numbers: 73.23.-b, 73.40.Gk, 73.50.Td 



One-dimensional (ID) interacting fermionic systems 
are described as Luttinger liquids (LL) whose exper- 
imental manifestations include carbon nanotubes, semi- 
conductor, metallic and polymer nanowires, as well as 
quantum Hall edges. While equilibrium LL have been 
extensively explored, there is currently a growing inter- 
est in non- equilibrium phenomena on the nanoscale. The 
major obstacle for a theoretical study of quantum wires 
out of equilibrium is that bosonization, a technique most 
suitable for analytic treatment of LL, has been so far 
restricted to equilibrium situation. 

In this work we develop a bosonic theory for inter- 
acting ID fermions under generic non-equilibrium con- 
ditions. The configuration we have in mind is shown in 
Figdl electrons with distributions n,,(e) {rj = R,L labels 
right- and left-movers) are injected into the LL wire from 
two non- interacting electrodes. A particularly interesting 
case involves double-step distributions [5| . Our goal is to 
calculate the electron Green functions (GFs) G<(r) mea- 
surable in tunneling spectroscopy experiments [3, IJ] . 

Our microscopic description of the problem begins with 
Keldysh action in terms of fermionic fields ip(t,x) [5j], 

Sec['lp]^-Y^ dt dxg{x){pr,P-,T+ PrjPn) ■ (1) 
^ Jc J 

Here pri = ^lii^ri are density fields, dri = dt + r/udx, v is 
the Fermi velocity, g{x) is a spatially dependent interac- 
tion strength, and 77 = ±1 stands for right/left moving 
electrons. The essence of bosonization is a reformulation 
of the theory in terms of bosonic (density) fields. The 
interacting part of the action See is already expressed in 
terms of density modes pri- The challenge is to bosonize 
the free part of the action, where information concern- 
ing the state of the non-interacting fermionic system is 



FIG. 1: Schematic view of a LL conductor connected to leads 
with two different incoming fermionic distributions. The LL 
interaction parameter K{x) is also shown; the dashed line 
corresponds to its sharp variation at the boundaries. 

encoded. The bosonic counterpart of the free part of the 
action is found by requiring that it reproduces correlation 
functions of non-equilibrium free Fermi gas Q, 



(2) 



Here we have performed a rotation in Keldysh space, de- 
composing fields into classical and quantum components, 
p, p = {p+ ± p-)/y/2, where the indices -I- and — refer to 
the two branches of the Keldysh contour; HJJ is the ad- 
vanced component of polarization operator, and Zrilxrj] 
is a partition function of free fermions moving in the field 



Xn = 



H° 



Pn- 



(3) 



Expansion of Z^ilxri] in Xvi^^^) generates an infi- 
nite series, zlnZ^[xr,] = I]n(-l)"^^X»7"'5';"/'^7 governed 
by irreducible fermionic density correlation functions, 
Sr,niti,xi; . . .■,tn,Xn) = {(p„iP2 ■ ■ ■ P^n)) , representing 
cumulants of quantum noise 8] . At equilibrium the ran- 
dom phase approximation (RPA) is exact 0|, iS^„ = 
for 71 > 2, and one obtains the Gaussian theory of con- 
ventional bosonization. In a generic non-equilibrium sit- 
uation, the bosonic action contains terms of all orders 
with no small parameter; the idea to proceed analyti- 
cally in a controlled manner may seem hopeless. This, 



however, is not the case: non-equilibrium bosonization is 
an efficient framework in which the functional integration 
can be performed exactly. Two observations are crucial 
here. First, depends only on the quantum compo- 
nent p, so that the action, Eq.®, is hnear with respect 
to the classical component p of the density field. Hence 
the integration with respect to p can be performed ex- 
actly, generating an equation that fixes p. As a result, 
the many-body field theory is mapped onto a problem 
which has a level of complexity of quantum mechanics. 
Second, all vertex functions Sr^n determining Z^^p^ are 
restricted to the mass shell (uj — r\vq) with respect to all 
of their arguments. This will render the effective "quan- 
tum mechanics" one-dimensional. 

We apply this formalism first to free fermions 
away from equilibrium, and obtain the single-particle 
GFs G<^^{x,t) = z(V'J(0,0)V^(x,t)), Gl^(x,t) = 
— t)i/)t (0, 0)) by means of the bosonic theory. 
Since Go,r; depends otlt — t — rjxjv only, we set x = 
and < = r in the arguments of Go.r;. The fermionic op- 
erators are expressed in terms of the bosonic field in the 
standard way, ~ (A/27rw)^/2e''*''^^e''^''(^\ where 

Priix) = rjdx4>ij'2.ii and A is an ultraviolet cut-ofF. Af- 
ter integration over the classical component we find that 
quantum component p satisfies a continuity equation 



dtPr) + VvdxPr, = jit,x) 



(4) 



with the source term j = d{x)[S{t — r) — 6{t)]/^/2. The 
fermionic GF is obtained as [with Xn given by Eq.iP])] 



Go<„(t) 
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(5) 



The mass-shell nature of S'^„ implies that [xjj] depends 
only on the world-line integral 

/oo 
dixr,{i~t,T]vi). (6) 
-oo 

More specifically, we find 

Z^x^] = det[l + (e-^^"" - = A„[<5„(i)]. (7) 

For the free-fermion problem the phase Srj{t) = Xujr{t,0) 
where Wr{t,t) = 9{t — t) — 9{t — t — t) is a "window 
function" and A — 27r. Thus, Zfj[xri] — A,,T-(27r), where 
Ajjr (A) is the determinant ([7]) for a rectangular pulse [7| . 

Since firi is diagonal in energy space, while Sr^ is di- 
agonal in time space, they do not commute, making 
the determinant non-trivial. Determinants of the type 
([7]) have appeared in a theory of counting statistics [8|. 
Specifically, the generating function of current fluctua- 
tions k(A) = X^S^-oo e™'^Pn (where p„ is the probabihty 
of n electrons being transferred through the system in a 
given time window r) has the same structure as A.,,r(A). 
Taylor expansion of In k(A) around A = defines cumu- 
lants of current fluctuations. 
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FIG. 2: Analytic Anri^) vs. periodic k(A) continuation of 
the functional determinant. The value of A,,r at A = 27r 
determines the free-electron GF, while In K(27r) = in view 
of charge quantization. As an example, the equilibrium case 
is shown. Inset: contour of integration for the quasiclassical 
limit, Eq. ((Sj, of A^t(A) is deformed, since a singularity of 
the integrand crosses the real axis at A = vr. 

According to its definition, k(A) is 27r-periodic, which 
is a manifestation of charge quantization. Thus, k{2tt) = 
1 is trivial. On the other hand, we have found that the 
free electron GF is determined by the value of the func- 
tional determinant exactly at A = 27r. A resolution of this 
apparent paradox is as follows: the determinant A^^(A) 
should be understood as an analytic function of A in- 
creasing from to 27T. On the other hand, k(A) is non- 
analytic at the branching points A — ±7r,±37r, .... To 
demonstrate this, it is instructive to consider the equi- 
librium case (temperature T). Then the expansion of 
In A,,,- (A) in A is restricted to the A^ term (since RPA is 
exact). It is easy to check that the A = 27r point on this 
parabolic dependence correctly reproduces the fermion 
GF via Eqs. ([7]). As to the counting statistics In k(A), 
it is quadratic only in the interval [— tt, tt] and is period- 
ically continued beyond this interval, see Fig. [51 

The difference in analytical properties of k{X) and 
A^T-(A) becomes especially transparent if one studies the 
semiclassical (long-r) limit. 

In A„,[A] f de\n[l + (e"^^ - l)n,(e)] . (8) 

For small A the singularity of the integrand closest to 
the real axis is located near e = iirT. As A increases, 
the singularity moves towards the real axis, crosses it at 
X ~ TT and finally approaches e = — ittT as A ^ 27r (see 
inset of Fig[21). The integral for InK(A) is taken along 
the real axis, resulting in non-analyticity at A = tt and in 
zero value at A = 27r. On the other hand, the contour for 
In ArjT (A) is deformed to preserve analyticity, resulting in 
the long-T asymptotics InA^,- = —tTX'^ /Air. 

The value of A = 27r appearing in the bosonic represen- 
tation of the free-fermion GF Go,,, (r) has a clear physical 
meaning: a fermion is a 27r-soliton in the bosonic for- 



malism. A natural question to ask is whether values of 
A^T-(A) away from A = 27r are physically important. To 
see that this is indeed the case, consider the Fermi edge 
singularity (FES) problem. In this problem, an electron 
excited into the conduction band, leaves behind a local- 
ized hole, resulting in an s-wave scattering phase shift, 
Sq, of the conducting electrons In the mesoscopic 

context, the FES manifests itself in resonant tunneling 
experiments The problem is effectively described 



by chiral ID electrons interacting with a core hole that 
is instantly "switched on" . As was shown in , tak- 
ing into account the core hole in the bosonization ap- 
proach amounts to a replacement of e*"^ by e*^^~''°/'^-''^ 
in the boson representation of the fermionic operator. 
Within our non-equilibrium formalism, this implies a re- 
placement j — > (1 — 5olT^)j in Eq. Q. Performing the 
derivation as in the free-fermion case, we thus obtain the 
non-equilibrium FES GF for electrons with an arbitrary 
distribution n(e). 
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G^{t) = =FjAA^[2(7r - (5o)]/27rw(l ± ikr) 



(9) 



It is easy to check that at equilibrium this reproduces 
known results 13, 12|. For a non-equilibrium double-step 
distribution the FES problem was studied in ^IS*! within a 
complementary (fermionic) approach (separating of G{t) 
into a product of open line and closed loop contribu- 
tions) bearing an analogy with functional bosonization 
[l^ . While both approaches yield equivalent results, this 
equivalence is highly non-trivial 6]. Functional deter- 
minants of the type At-(A) can be efficiently evaluated 
numerically [l5| . 

We are now prepared to address the problem formu- 
lated in the beginning of the paper: an interacting quan- 
tum wire out of equilibrium, FiglT] We consider G'^(r) 
for the tunneling point (x = 0) located inside the in- 
teracting part; generalization to tunneling into a non- 
interacting region is straightforward. Integrating over 
the classical components as explained above, we ob- 
tain equations satisfied by the quantum components p,, 
of the density fields. 



dtPR + dx iv + —)pR 

ZTT 



2tt 
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ZTT ZTT 



= J 

= 0. 



(10) 



The solution of Eq. pO)) determines the phases 5n{t) ac- 
cording to Eqs. (O, which can be cast in the form 



5r,{t) = -27r\/2 77 Jim 

t — > — oc 



»)f (t+t) 



dx Priix, t) . (11) 



Remarkably, Eq. (fTTj) expresses the phase 5,j{t) affected 
by the electron-electron interaction, through the asymp- 
totic behavior of p{x,t) in the non- interacting parts of 
the wire (regions I and III in Fig[T|). The phases 5jj{t) 



FIG. 3: Phases Sr^{t) entering Eq. ITU) for the GFs for sharp 
[a),b); r — {1 — K)/{1 + K)] and adiabatic [c),d)] boundaries. 



determine the GFs via 



G|(r) = - 



An[SR{t)]AL[5L{t)] 



(12) 



where 7 = (1 - Kf /2K and K = [l + gjixvY^I'^ is the 
standard LL parameter in the interacting region. 

To explicitly evaluate (5,, it) for the structure of Fig. [U 
it is convenient to rewrite Eqs. (jlOp as a second-order 
differential equation for the current J — v{pj^ — Pl)j 



{lo^ + d^u^ix)d.^)J{Lj,x) ^0, x^O, 



(13) 



where u{x) — v{l + g{x)/'Kv)^^'^ is a spatially depen- 
dent plasmon velocity. Reflection and transmission of 
plasmons on both boundaries is characterized by the co- 
efficients r,,, t,i (r^ + ~ 1); for simplicity, we assume 
them to be constant over a characteristic frequency range 
uj ~ T^^. Subsequently Sjj{t) is found to be a superpo- 
sition of rectangular pulses, S^iit) — J2'^^o^v,n'^Ti^'^n)i 
where t„ (n+ 1/2 - l/2K)L/u, and 



"'-"'+\1-7]K)/Vk . (14) 



The phases Sjf{t) are shown in Fig. [3] for two limits of 
adiabatic (r^ = 0) and sharp [r,, = (1 — K)/{1 + K)] 
boundaries. 

In physical terms, Sri{t) characterizes phase fluctua- 
tions in the leads that arrive at the time interval [0, r] at 
the measurement point x — 0, governing the dephasing 
and the energy distribution of electrons encoded in the 
GFs G^{t). Up to inversion of time, one can think of 
5r,{t) as describing the fractionalization of a phase pulse 



(electron-hole pair) injected into the wire at point x dur- 
ing the time interval [0,t]. This is closely related to 
the physics of charge fractionalization discussed earlier 
lil 17 1 . At the first step, the pulse splits into two with 
relative amplitudes (1 -I- K)/2 and (1 — K)J2 carried by 
plasmons in opposite directions, cf. Ref. [l6|. As each 
of these pulses reaches the corresponding boundary, an- 
other fractionalization process takes place: a part of it is 
transmitted into a lead, while the rest is reflected. The 
reflected pulse reaches the other boundary, is again frac- 
tionalized there, etc. Let us stress an important differ- 
ence between boundary fractionalization of transmitted 
charge [l^ and that of dipole pulses. While in the former 
case the boundaries can always be thought as sharp, in 
the present problem the way K{x) is turned on is cru- 
cially important for reflection coefhcients at a; ~ t^^. 

For T L/u the coherence of plasmon scattering may 
be neglected and the result splits into a product 



(15) 



with each factor representing a contribution of a single 
phase pulse (5,,,n(i) — (5,,,„w,-(t, 0). For the "partial equi- 
librium" state (when nR{t) and 71^(1) arc of Fcrmi-Dirac 
form but with different temperatures) the functional de- 
terminants are gaussian functions of phases, reproducing 
earlier results of functional bosonization ■ 

We now apply our general result (fT5| to the "full non- 
equilibrium" case, when n^(e) have a double step form, 
nr)(e) = ajj9{Vjj — e) -I- (1 — a^)0(— e). The split zero-bias- 
anomaly dips are then broadened by non-equilibrium de- 
phasing rate, 1/r^ @jl3l- The latter can be found from 
the long-time asymptotics of G'^(r), employing Eq. ([8]). 
We obtain l/r^ = + where 1/r^''' is the 

contribution to dephasing of 77 fermions governed by the 
distribution of rj' fermions. Focussing on the adiabatic 
limit, we find these dephasing rates to be 



-£^ln 
27r 



4a,j(l — a,j) sin^ 



7r(l + T]K 



2VK 



(16) 

see Fig. 01 Two remarkable features of this result should 
be pointed out. First, for certain values of the inter- 
action parameter K (different for 77 = R,L) the de- 
phasing rates ^/t^^ vanish. Second, consider the weak- 
interaction regime, 7^1. We then obtain ~ 
TT-yeVLaLil - ul) and 1/t^« ~ 71(778)6^^0^(1 - gr). 
This should be contrasted with RPA which predicts equal 
l/r^^ and l/r^^ [2,]. While l/r^^ agrees with the RPA 
result, is parametrically smaller (suppressed by 

an extra factor of 7) . The reason for this failure of RPA 
is clear from our analysis. For a weak interaction the 
contributions of R and L movers to Gr are given by the 
functional determinants A^t-(^^) with phases (for an adi- 
abatic barrier) 5^ — (1 — K)tt and Sr ~ Tr{l + K). While 
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FIG. 4: Dephasing rates and t^'" as function of LL pa- 
rameter K for the adiabatic case and double step distributions 
with an ~ aL = 1/3. 

the contribution of the small phase 5l is captured cor- 
rectly by RPA, a small-i5 expansion of \n. /S.rt-{5r) fails 
for large 5r (apart from equilibrium and "partial equi- 
librium" where In A^t-(<5) cx (S^.) 

To conclude, we have developed a non-equilibrium 
bosonization technique and used it to find GFs for tun- 
neling spectroscopy of a LL conductor. Our solution is 
presented in terms of functional determinants of single 
particle operators. It elucidates the relation between the 
charge fractionalization on one hand, and energy redistri- 
bution and non-equilibrium dephasing on the other hand. 
We have calculated dephasing broadening of the ZBA 
dips for double-step distributions. Our approach has var- 
ious further applications, including spinful electrons, in- 
terferometry, noise, and many-body entanglement [6[. 
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